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Abstract 

This paper shows that the left-invariant geodesic flow on the symplectic 
group relative to the Frobenius metric is an integrable system that is not con¬ 
tained in the Mishchenko-Fomenko class of rigid body metrics. This system 
may be expressed as a flow on symmetric matrices and is bi-Hamiltonian. This 
analysis is extended to cover flows on symmetric matrices when an isomorphism 
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with the symplectic Lie algebra does not hold. The two Poisson structures as¬ 
sociated with this system, including an analysis of its Casimirs, are completely 
analyzed. Since the system integrals are not generated by its Casimirs it is 
shown that the nature of integrability is fundamentally different from that ex¬ 
hibited in the Mischenko-Fomenko setting. 


1 Introduction 

This paper continues the analysis, begun in Bloch and Iserles [2006], of the set of 
ordinary differential equations 


X = [X^,N], (1.1) 

where X £ Sym(n), the linear space of ri x n symmetric matrices, X denotes the 
time derivative, N £ 5o(n), the space of skew symmetric n x n matrices, is given, 
and where initial conditions 2f(0) = Xq £ Sym(n) are also given. 

It is easy to check that [X^,A^] £ Sym(n), so that if the initial condition is in 
Sym(n), then X{t) £ Sym(n) for all t. Also, because of the straightforward identity 
[A^, a] = [X, XN + NX], this equation may be rewritten in the Lax form 

X = [X,XN + NX], (1.2) 

again with initial conditions A(0) = Xq £ Sym(n). 

We show below that this system may be viewed as a Lie-Poisson system on the 
dual of the symplectic Lie algebra if N invertible, in which case it is geodesic, and on 
the dual of a more general Lie algebra on symmetric matrices for arbitrary N. The 
system is bi-Hamiltonian and is not in the Mischenko-Fomenko class of integrable 
(geodesic) rigid body systems (Mishchenko and Fomenko [1976]). Despite this, we 
prove that it is integrable on the generic symplectic leaf of the corresponding phase 
space if N is invertible or of nullity one. We use the Lax pair with parameter 
found in Bloch and Iserles [2006] to find a class of integrals that we show are in 
involution using the bi-Hamiltonian structure and the technique of Morosi and Piz- 
zocchero [1996]. Independence is proved directly since the method in Mishchenko 
and Fomenko [1976] does not apply to this system, even though the integrals are 
obtained from the Casimirs with a shifted argument. Indeed, this system appears 
to fundamentally different from completely integrable systems either of rigid body 
or Toda type (on symmetric matrices). 

If N is not invertible, there is no isomorphism of the Lie algebra induced by 
N with the symplectic Lie algebra. We extend our analysis of the system to this 
case and study the Poisson geometry of the dual of this Lie algebra determining 
the generic leaves and the Casimir functions of both Poisson structures relative to 
which the system (1.1) is bi-Hamiltonian. 

We want to emphasize that the system (1.1) (or (1.2)) for N invertible is thus a 
new integrable geodesic flow of a left invariant metric on the Lie group Sp(n,M). So 
far the only known left invariant metrics whose geodesic flows are integrable on the 
Lie group Sp(n,M) are the rigid body metrics of Mishchenko and Fomenko [1976]. 
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Finding integrable geodesic flows on Lie groups for left invariant metrics that are 
not of rigid body type is a daunting task. System (1.1) is the only one known to 
us on any semisimple Lie algebra with the exception of 5o(4), where we review the 
stituation below. 

Even for the case of SO(4) there is only one known geodesic flow that is not of 
rigid body type (see Mishchenko [1970]; Manakov [1976]; Mishchenko and Fomenko 
[1976]; Ratiu [1980] for the definition of such metrics). There are three integrable 
cases of left invariant metrics for geodesic flow on SO(4): the metric used in Manakov 
[1976] (which goes to the Clebsch case by contraction to the Euclidean group) and 
two other cases that correspond to left invariant metrics that are not diagonal in the 
standard basis of 5o(4). The first one is obtained from deformation of the classical 
Lyapunov-Steklov integrable case on SE(3) by deforming the Lie algebra 5e(3) to 
so (4); the integrability of the corresponding system is due to Borisov, Mamaev, and 
Sokolov [2001]. The last case has a fourth quartic constant of the motion and is the 
genuinely new integrable geodesic case found by Adler and van Moerbeke [1986]; 
a Q 2 Lax pair for this system was given in Reyman and Semenov-Tian-Shansky 
[1986]. Sokolov [2001] showed that these two cases are not linearly equivalent. 
The rigid body metric used by Manakov [1976] is the only algebraic completely 
integrable case for a left invariant metric that is diagonal in the standard basis of 
50 (4) (Adler and van Moerbeke [1982]; Haine [1984]). The state of the art regarding 
these systems is contained in Theorem 8.3, page 270, of Adler, van Moerbeke, and 
Vanhaecke [2004]: in a certain large class of metrics (non-degenerate half-diagonal 
metrics with some weight homogeneity conditions) these three cases are the only 
algebraically completely integrable geodesic flows. Whether these three cases are 
the only algebraically completely integrable geodesic flows in the class of all left 
invariant metrics is still an open question. See Sokolov [2002] for a review and 
references of what is known about a related system, the Kirchhoff case of the motion 
of a rigid body in an ideal fluid. 

The structure of the paper is as follows. In Section 2 we consider the Lie algebra 
structure on symmetric matrices induced by N and the special case of the isomor¬ 
phism to sp(M”'). In Section 3 we analyze the bi-Hamiltonian structure as well as 
the symplectic leaves and Casmirs of both structures. In Section 4 we compare our 
system with the sectional operator systems of Mischenko and Eomenko and conclude 
that (1.1) is not in this family, thereby showing that it is a new geodesic flow that 
is not of rigid body type on the Lie group Sp(nM). In Section 5 we analyze the Lax 
pair with parameter and find a family with the right number of integrals of motion 
that is a candidate for Liouville integrability. In Section 6 we prove involution of the 
integrals using the bi-Hamiltonian structure. In Section 7 we analyze indepedence 
and finally we discuss some future work in Section 8. 

2 The Lie Algebra and the Euler-Poincare Form 

We can regard as a Poisson tensor on M” by defining the bracket of two functions 
f,9 as 




( 2 . 1 ) 


2 The Lie Algebra and the Euler-Poincare Form 


4 


The Hamiltonian vector field associated with a function h (with the convention 
that f{z) = Xh{z) ■ Vf{z) = {/, h} {z)) is given by 


Xhiz) = NXhiz) 


( 2 . 2 ) 


as is easily checked. 

For each X G Sym(n) define the quadratic Hamiltonian Qx by 

Qx{z) := -z"^Xz, z G M”. 

Let Q := {Qx | X G Sym(n)} be the vector space of all such functions. Note that 
the map Q : X e Sym(n) i— Qx G Q is an isomorphism. 

Using (2.2) it follows that the Hamiltonian vector field of Qx has the form 

Xq^{z) = NXz. (2.3) 

Next, we compute the Poisson bracket of two such quadratic functions. 

Lemma 2.1. For X,Y G Sym(n), we have 

{Qx, Qy}n = Q[x,y]n^ (2-4) 

where [X,Tjjv = XNY — YNX G Sym(n). In addition, Sym(n) is a Lie algebra 
relative to the Lie bracket [-j-Jat. Therefore, Q : X € (Sym(n), [•,-jAr) i—> Qx G 
® Lie algebra isomorphism. 

Proof. Using (2.1) we have 

{Qx, Qy}n (^) = (VQx) {zfN (VQy) (z) = (Xzf NYz = z^XNYz 
= {XNY - YNX) z = Q[x,y]A^)- 

Recall that the notation Qy is reserved only for symmetric matrices V. Since 
X,Y G Sym(re) implies that [X,Y]x = XNY — YNX G Sym(re) we can write 
Q[x,Y]j.f in the preceding equation. 

The bracket [•, -Jat on Sym(n) is clearly bilinear and antisymmetric. The Jacobi 
identity is a straightforward direct verification. ■ 

It is a general fact that Hamiltonian vector fields and Poisson brackets are related 

by 

= (2.5) 

where the bracket on the left hand side is the Jacobi-Lie bracket. Thus, it is natural 
to look at the corresponding algebra of Hamiltonian vector fields on the Poisson 
manifold (M”, {•, -jAf) associated to quadratic Hamiltonians. If we take / = Qx and 

g = Qy, with Xf = NX and Xg = NY, and recall that the Jacobi-Lie bracket 

of linear vector fields is the negative of the commutator of the associated matrices, 
then we have the following result. 
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Proposition 2.2. Equations (2-4) and (2.5) imply 

N[X,Y]n = [NX,NY]. (2.6) 

This can, of course, be easily verified by hand. 

Letting C7i denote the Lie algebra of linear Hamiltonian vector fields on M"' 
relative to the commutator bracket of matrices, (2.6) states that the map 

A G (Sym(n), [•, -J^) ^ NX ^ {CH, [•, •]) 

is a homomorphism of Lie algebras^. 

Invertible Case. If N is invertible, then this homomorphism is an isomorphism. 
In addition, the non-degeneracy of N implies that n is even and that is a sym- 
plectic vector space relative to the symplectic form defined by N~^. Therefore, the 
Lie algebra {CTL, [•, ■]) is isomorphic to the Lie algebra sp(M”, A“^) of linear sym¬ 
plectic maps of M"' relative to the symplectic form that is, to the classical Lie 
algebra 5p(n,M). Note that this means that {NX)"^N~^ + N~^{NX) = 0. 

We summarize these considerations in the following statement. 

Proposition 2.3. Let N G so(n). The map Q : X € (Sym(n), [•,-Jat) i—> Qx G 
(Qji'rlAf) a Lie algebra isomorphism. The map X G (Sym(n), [-j-Iat) e-> NX G 
{CTL, [•, •]) is a Lie algebra homomorphism and if N is invertible it induces an iso¬ 
morphism o/(Sym(n), [•,-jAr) with 5p{n,'R). 

The Euler-Poincare Form The Euler-Poincare form for the equations can be 
derived as follows. Identify Sym(re) with its dual using the the positive definite inner 
product 

{{X, Y)) := trace {XY ), for A, T G Sym(n). (2.7) 

Remark. The inner product ((A, T)) is not ad-invariant relative to the A-bracket, 
but another one, namely kx{X,Y) := trace(AAAy) is invariant, as is easy to 
check. 

Define the Lagrangian I : Sym(n) —M on the Lie algebra (Sym(n), [•, -jAf) by 

1{X) = i trace (A^) = i trace (AA^) =: ^ ((A, A)) . (2.8) 

Proposition 2.4. The equations 

X = [X‘^,N] (2.9) 

are the Euler-Poincare equations^ corresponding to the Lagrangian (2.8) on the Lie 
algebra (Sym(n), [■,-]n)- 

^We thank Gopal Prasad for suggesting isomorphisms of this type; they are closely related to 
well-known properties of linear Hamiltonian vector fields, as in Marsden and Ratiu [1994], Propo¬ 
sition 2.7.8. 

^For a general discussion of the Euler-Poincare equations, see, for instance, Marsden and Ratiu 
[1994]. 
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Proof. Recall that the general (left) Euler-Poincare equations on a Lie algebra q 
associated with a Lagrangian ^ : g —> M are given by 

j^Dm = adiDm, 

where Dl{^) G g* is the Frechet derivative of I at Equivalently, for each fixed 
?7 G g, we have 

±Dm-rj = Dm-[C,r]]. ( 2 . 10 ) 

In our case, letting ^ = X and r] = Y arbitrary, time-independent, equations (2.10) 
become 

j^{{X,Y)) = {{X,[X,Y]n)) 

= {{X, XNY - YNX )); 

that is, 

trace = trace {X{XNY - YNX)) 

= trace {{X^N - NX‘^)Y) , 

which gives the result. ■ 

General Case—Noninvertible N. We next determine the structure of the Lie 
algebra (Sym(n) A;-]n) for a general skew-symmetric matrix N. The point of 
departure is the fact that if N is nondegenerate, then X G (Sym(n), [•,'Jat) i--> 
NX G {CH, [■,■]) = (5p(]R"', X“^), [•, ■]) is a Lie algebra isomorphism. Recall that 
if M” has an inner product, which we shall take in what follows to be the usual 
dot product associated to the basis in which the skew-symmetrix matrix N is 
given, and L : M” —> M” is a linear map, then M"" decomposes orthogonally as 
M"" = imL^ 0 ker L. Taking L = N in this statement and recalling that N"’" = —N, 
we get the orthogonal decomposition M"' = imX © kerX. Let 2p = rankX and 
d := n — 2p. Then N := X|ijn 7 v : imX —> imX defines a nondegenerate skew 
symmetric bilinear form and, by the previous proposition, (Sym(2p), [■,-J^y) is iso¬ 
morphic as a Lie algebra to (5p(M^^, iV“^), [■, •]). In this direct sum decomposition 
of M"', the skew- symmetric matrix N takes the form 


where iV is a {2p) x (2p) skew-symmetric nondegenerate matrix. 

The Lie algebra (Sym(2p), [■, -J^y) acts on the vector space X\[ 2 p)xd of {2p) x d 
matrices (which we can think of as linear maps of ker X to imX) by S' • A := SNA, 
where S G (Sym(2p), [• -J^y) and A G ■M( 2 p)xd- Indeed, if S,S' G Sym(2p) and 
A G XI( 2 p)xd, then 

[S, S'jjy ■ A = (SNS' - S'NS)NA = SNS'NA - S'NSNA 
= S ■ (S' ■ A) - S' ■ (S ■ A). 


( 2 . 11 ) 
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Now form the semidirect product Sym(2p) (D A^( 2 p)xd- Its bracket is defined by 

[{S, A), (5', A')] = ([5, S%, S-A'-S'-A) 

= {SMS' - S'NS, SNA' - S'NA) (2A2) 

for any S, S' G Sym(2p) and A, A' G ■M( 2 p)xd- 

Next, define the Sym(d)-valued Lie algebra two cocycle 

C : Sym{2p)®M^2p)xd x Sym{2p)®M^2p)xd ^ Sym(d) 

by 

CHS, A), {S', A')) := A'^NA' - {A'fNA (2A3) 

for any S, S' G Sym(2p) and A, A' G A4(^2p)xd- The cocycle identity 

C([(5, A), (S', A')], {S", A")) + Ci[{S', A'), {S",A")], {S, A)) 

+ C([(5",A"),(5,A)],(5',A')) = 0 


for any S,S',S" G Sym(2p) and A, A', A" G Ni(2p)xd is a straightforward verifica¬ 
tion. Now extend Sym(2p) (s)Ad( 2 p)xd by this cocycle. That is, form the vector 
space (Sym(2p) (s) A4(2p)xd) ® Sym((i) and endow it with the bracket 

[{S, A, B), (S', A', B')f := (SNS' - S'NS, SNA' - S'NA, A^NA' - (A'fNA) 

(2.14) 

for any S, S' G Sym(2p), A, A' G Ni{ 2 p)xdi and B,B' € Sym(d). 

Proposition 2.5. The map 

T : ((Sym(2p) (s) A4(2p)xd) © Sym(d), [•, ^ (Sym(n), [■, -(n) 


given by 


^{S,A,B) 


is a Lie algebra isomorphism. 


■ S A 
A^ B 


(2.15) 


Proof. It is obvious that T is a vector space isomorphism so only the Lie alge¬ 
bra homomorphism condition needs to be verified. So, let (S, A, B), (S', A', B') G 
(Sym(2p) (s) A4(2p)xd) © Sym((i) and compute 


A>{[{S, A, B), (S', A', B')]) = Ai{SNS' - S'NS, SNA' - S'NA, A^NA' - (A'fNA) 

SNS'-S'NS SNA'-S'NA_ 

(SNA' - S'NA)'^ A^NA' - (A')'^NA 



'N O' 

S' A'' 


S' A'' 

'N O' 

' s 

A 


0 0 

(A')'^ B' 


(A')'^ B' 

0 0 

A^ 

B 


= [A>(S,A,B),A^(S',A',B')U 


as required. 
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3 Poisson Structures 


Identifying Sym(n) with its dual using the inner product (2.7) endows Sym(n) with 
the the (left, or minus) Lie-Poisson bracket 


{f,9}N (^) = - trace 


x(vf{X)NVg{X) - XgiX)NXf{X)^], 


(3.1) 


where V/ is the gradient of / relative to the inner product ((•,•)) on Sym(n). It 
is easy to check that the equations X = are Hamiltonian relative to the 

function I defined in (2.8) and the Lie-Poisson bracket (3.1). 


Later on we shall also need the frozen Poisson bracket 

if, gj^j, {X) = - trace (vfiX)NVg{X) - Xg{X)NVf{X)) . (3.2) 

It is a general fact that the Poisson structures (3.1) and (3.2) are compatible in the 
sense that their sum is a Poisson structure (see e.g. Exercise 10.1-5 in Marsden and 
Ratiu [1994]). 

For what follows it is important to compute the Poisson tensors corresponding 
to the above Poisson brackets. Recall that the Poisson tensor can be viewed as a 
vector bundle morphism B : r*(Sym(n)) —> r(Sym(n)) covering the identity. It is 
defined by B{dh) = {■,h}]\f for any locally defined smooth function h on Sym(n). 
Since Sym(n) is a vector space, these bundles are trivial and hence the value Bx at 
X G Sym(n) of the Poisson tensor R is a linear map Bx ■ Sym(n) —> Sym(re) by 
identifying Sym(n) with its dual using the inner product ((•, •)). 

Proposition 3.1. Denote the value at X £ Sym(n) of the Poisson tensors corre¬ 
sponding to the Lie-Poisson (3.1) and frozen (3.2) brackets by Bx and Cx, respec¬ 
tively. Then for any Y £ Sym(n) we have 

BxiY) = XYN - NYX (3.3) 

Cx{Y) = YN - NY. (3.4) 

Proof. Let / and g be locally defined smooth functions on Sym(re). The definition 
of Bx gives 


{{Vf{X),Bx{Ng{X))) 


= {f,g}NiX) 


= — trace 


X 


(xf{X)NXg{X) - Xg{X)NVf{X)) 
Xf{X)(^XVg{X)N - NXg{X)X^ 
{{XfiX),XXg{X)N - NXg{X)X)), 


= trace 


which implies (3.3) since any Y £ Sym(n) is of the form Xg{X), where g{X) 








3 Poisson Structures 


9 


{{X,Y)). Similarly, the definition of Cx gives 


{{Vf{X),Cx{V9{X))) 


= {f,9}FN{X) 

= - trace (xf{X)NXg{X) - Xg{X)NXf{X)) 

= trace Vf(X) (Xg{X)N - NXg{X)'^ 

= {{Vf{X),Xg{X)N - NXg{X))), 


which proves (3.4). ■ 

Proposition 3.2. Let n = 2p + d, where 2p = rankiV. The generic leaves of the 
Lie-Poisson bracket {-j-Iat are 2p{p + d)-dimensional. 


Proof. As in the proof of Proposition 2.5, we orthogonally decompose M” = im A^0 
ker N so that iV = A^| im : im > im A^ is an isomorphism. In this decomposition 
the matrix N takes the form 


and, according to the isomorphism ^ in Proposition 2.5, the matrix X can be written 
as 


A = 


S 


A 

B ’ 


where S € Sym(2p), B € Sym(d), and A € A4(2p)xd- Therefore, if 


U C 


D 


G Sym(n) 


with U G Sym(2p), D G Sym.{d), C G M.[ 2 p)xdi the Poisson tensor of the Lie-Poisson 
bracket {•, takes the form (see Proposition 3.1) 

Bx{Y) = XYN - NYX 



r u 

cl 

'N O' 


'N O' 

r u 

cl 

\ S 

Al 



D 

0 0 


0 0 


D 

A^ 

B 


_ 'SUN - NUS + AC^N - NCA^ -NUA-NCB 
A^UN + BC'^N 0 

Since N is invertible, the kernel of Bx '■ Sym(n) —> Sym(n) is therefore given by all 
U G Sym(2p), C G Sym(d), and C G Xi( 2 p)xd such that 

SUN - NUS + AC'^N - NCA^ = 0 and UA + CB = 0. 


To compute the dimension of the maximal symplectic leaves, we assume that the 
matrix X is generic. So, supposing that B is invertible, we have C = —UAB~^ and 

{S - AB-^A^) UN -NU{S- AB-^A^) = 0. 
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Since S—AB~^A^ E Sym(2p) is given, this condition is identical to the vanishing 
of the Poisson tensor on the dual of the Lie algebra (Sym(2p), [■, -j^) evaluated at 
S — But N is invertible so, according to Proposition 2.3, this Lie algebra 

is isomorphic to 5p(2p, M) whose rank is p. Therefore, the kernel of the map 

U E Sym(2p) ^ (5 - AB-^A^) UN - NU {S - AB-^A^) E Sym(2p) 

for generic S — AB~^A^ has dimension p. 

Since C = —UAB~^ is uniquely determined and D E Sym((i) is arbitrary, we see 
that the dimension of the kernel of Bx for generic X has dimension p + d{d + l)/2. 
Thus the dimension of the generic leaf of the Lie-Poisson bracket {•, -Iat is 

^{2p + d){2p + d + 1) - p - ^d{d + 1) = 2p{p + d) 

as claimed in the statement of the proposition. ■ 

Proposition 3.3. All leaves of the frozen Poisson bracket {•, ■}fn 

(i) 2p{p + d)-dimensional if N is generic, that is, all its non-zero eigenvalues are 
distinct, and 


(ii) p{p + 1 + 2d)-dimensional if all non-zero eigenvalue pairs of N are equal. 

Proof. Proceeding as in the proof of the previous proposition and using the same 
notation for N, X, and Y, the Poisson tensor of the frozen bracket takes the form 


Cx{Y) = YN -NY = 


■ u 

c 


'N O' 


'N O' 

■ u 

c 


D 


0 0 


0 0 


D 


'UN - NU NC 
C^N 0 


Thus, since N is invertible, the kernel of Cx is given by all U E Sym(2p), D E 
Sym(d), C E M.{2p)xd such that (7 = 0 and UN — NU = 0. 

Since N is non-degenerate, there is an orthogonal matrix Q such that 


N = Q^ 


0 

-V 


V 

0 


Q, 


where V = diag(ui,..., Vp) and u* E M, Uj / 0 for all i = 1, ... ,p. Therefore, 


0 = UN - NU = UQ 

0 V 


■ 0 V' 


' 0 

V' 

-V 0 

Q - Q 

-V 

0 


QU 


= QUQ 


0 V 


is equivalent to 


U 


-V 

0 

-V 

■ 0 

V' 


' 0 

V' 

-V 

0 


-V 

0 


QUQ^ Q 


U = 0 


(3.5) 


























3 Poisson Structures 


11 


where tj := QUQ^ G Sym( 2 p). Write 


- t/n C/i2 

[UL U22_ 

with Uii and U22 symmetric pxp matrices and U12 an arbitrary pxp matrix. Then 
( 3 . 5 ) is equivalent to 

U22 = VUuV-^ = V-^UiiV and = -V-^U^V = -VU^V-^. ( 3 . 6 ) 


(i) Assume now that Vi 7^ Vj if i 7^ j. Since VUnV ^ = V ^UnV is equivalent 
to V‘^UiiV~‘^ = C/11, it follows that 

for all z,j = l,...,p, 


where ttii,ij are the entries of the symmetric matrix C/n. Since the fraction on the 
left hand side is never equal to one for i ^ j, this relation implies that = 0 for 
all i 7 / j. Thus Uu is diagonal and U 22 = Uu. A similar argument shows that C /12 
is diagonal. However, then it follows that C /12 = ~Ui 2 which implies that C /12 = 0. 
Therefore, the kernel of the map U ^ UN — NU is p-dimensional. 

Concluding, the dimension of every leaf of the frozen Poisson structure equals 
\{2p + d){2p + d + 1) - p - \d{d + 1) = 2p{p + d). 

(ii) The other extreme case is when Vi = Vj =: v for all i,j = 1,... ,p. Then 
V = vl, where I is the identity matrix, and (3.6) becomes U 22 = Uu, C/)^ = 
—Ui 2 - Therefore, the kernel of the map U 1 —> UN — NU has dimension equal to 
ip(p + 1) + y{p-l) 

Concluding, the dimension of every leaf of the frozen Poisson structure equals 
^{2p + d){2p + d + 1) - p"^ - ^d{d + 1 ) = p{p +l + 2d). ■ 


Proposition 3.4 (Casimir Functions). Let the skew symmetric matrix N have 
rank 2p and size n := 2p + d. Choose an orthonormal basis in which N is 

written as 


N = 


0 

-V 

0 


V 0 
0 0 
0 0 


where V is a real diagonal matrix whose entries are vi,... ,Vp. 


(i) If Vi 7/ Vj for all i 7^ j, the p + d{d + l)/2 Casimir functions for the frozen 
Poisson structure ( 3 . 2 ) are given by 

Cp{X) = trace(£'jA), i = 1,... ,p + ^d{d + 1 ), 

where Ei is any of the matrices 


Skk 

0 

o' 


'0 

0 

0 ■ 

0 

Skk 

0 

5 

0 

0 

0 

_ 0 

0 

0 


0 

0 

Sab 
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Here S^k is the p x p matrix all of whose entries are zero except the diagonal 
(k, k) entry which is one and Sab is the d x d symmetric matrix having all 
entries equal to zero except for the (a, b) and {b, a) entries that are equal to 
one. 


(ii) If Vi = Vj for all i, j = 1,. .. ,p, the p^ + d{d + l)/2 Casimir functions for the 
frozen Poisson structure (3.2) are given by 

Cp{X) = trace(£'jX), i = 1,... ,p‘^ + + 1), 

where Ei is any of the matrices 


Ski 

0 

o' 


0 

Am 

o' 


'0 

0 

0 ' 

0 

Ski 

0 


-Aki 

0 

0 


0 

0 

0 

0 

0 

0 


0 

0 

0 


0 

0 

Sab 


Here Ski is the p x p symmetric matrix having all entries equal to zero except 
for the {k, 1) and {I, k) entries that are equal to one and Aki is the p x p skew 
symmteric matrix with all entries equal to zero except for the {k, 1) entry which 
is 1 and the {l,k) entry which is —1. 


(iii) Denote 


N = 


■ 0 

V' 


1 

o 

-V 

0 

and N = 

0 0 


The p + d{d + l)/2 Casimir functions for the Lie-Poisson bracket {•, - are 
given by 


C^{X) = ^ trace 

ZtKi 


(Xiv) 


2k 


, for k = l,...,p 


and 


C^(X) = trace(XSfc), for k = p + 1,... ,p-\ — d{d+l), 


where Ek is any matrix of the form 


0 0 0 
0 0 0 
0 0 Sab 


In the special case when N is full rank the Casimirs are just 


C’^iX) 


1 

— trace 
2k 




for = 


Proof. To prove (i), recall from Proposition 3.3(i) that the kernel of the Poisson 
tensor Cx has dimension p + ^d{d + 1). Moreover, if E belongs to this kernel, then 
the linear function given by X trace(£'X) has gradient E, which is annihilated 
by the Poisson tensor Cx- Thus all Cp are Casimir functions. Since the gradients 
of all these functions are the p + ^d{d + 1) matrices in the statement which are 
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obviously linearly independent, it follows that the functions Cp form a functionally 
independent set of Casimir functions for the frozen bracket {•, -Iftv- 
Part (ii) has an identical proof. 

In the proof of (iii) we do not need the detailed 3x3 block decomposition of 
N and X and shall use exclusively the 2x2 block decomposition, where the (1,1) 
block has size {2p) x (2p). Consider first the functions C^{X) for k = 1,... ,p. Note 
that VC^{X) = NXN ■ ■ ■ NXN (with {2k — 1) factors of X) and hence (3.3) gives 


(3.7) 


Bxi^C^iX)) = X NXN ■ ■ ■ NXN )N -N( NXN ■ ■ ■ NXN X. 


Note firstly that in the case N is invertible this is just XN ■ ■ ■ NX — XN ■ ■ ■ NX 
which is clearly 0. 

Now consider the general case. We first observe that 


NN = NN = 


I 0 
0 0 


(3.8) 


The product of the last four factors in the first term of equation (3.7) is thus 

NXNN = 


'N-^S O' 
0 0 


Similarly, the product of the first four factors of the second term of (3.7) is 

NNXN = 


SN-^ 

0 


Continuing the multiplication in both terms in this fashion (always taking a group of 
three consecutive factors from the right and left, respectively) we see that both terms 
have only nonzero (1,1) blocks which are identical and equal to SN~^SN~^....N~^S. 
Thus, again, (3.7) is identically zero. 

However, sp(2p, M) is identified with the subalgebra consisting of the (1,1) blocks 
of elements of Sym(n) (see Proposition 2.5). The isomorphism S € Sym(2p) i—> 
NS G 5p(2p, M) given in Proposition 2.3 identifies the basis of p Casimirs in the 
dual of 5p(2p, M) (given by the even traces of the powers of a matrix) with the 
functions S i—> trace [(S'iV“^)^^] /2k. Therefore the functions for k = 1,... ,p 
given in the statement of the proposition are functionally independent Casimirs for 
the Lie-Poisson bracket of Sym(n). 

To see that the remaining functions C^{X) = trace(XPl^) are Casimirs observe 
that in this case 

XC^{X) = ° 


0 Sab 


and 


Bx{VC\X)) = 


r 5 

Al 


'0 0 ■ 

'N O' 


'N O' 

'0 0 ■ 

\ ^ 

^1 


B 


0 Sab 

0 0 


0 0 

0 Sab 


B 


= 0 . 
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Since the matrices Sab span the symmetric k x k matrices, these Casimirs are func¬ 
tionally independent. The two sets of Casimirs are also independent taken together, 
since each set depends only on a subset of independent variables and these two sets 
of variables are disjoint. We have thus obtained p + d{d + l)/2 Casimirs, which is 
the codimension of the generic leaf thus proving that they generate the space of all 
Casimir functions of the Lie-Poisson bracket. ■ 


The equations in the degenerate case. If N is degenerate, representing it 
and the matrix X € Sym(n) as in Proposition 2.5, the equations X = [X'^,N] are 
equivalent to the system 

S=[S^ + A^A, N] 

A = -N{SA + AB) 

B = 0 


4 The Sectional Operator Equations 

This section shows that the flow (1.1) is not of the sectional operator type discussed 
in Mishchenko and Fomenko [1976]; in fact, this is the case already for 2x2 matrices 
with the canonical choice of N. 

Let 

0 1 
-1 0 


N = 


(4.1) 


and denote elements of Sym(2) by 

X = 


a b 
b d 


a,b,c ^ M. 


(4.2) 


One can readily check that a maximal Abelian subalgebra of Sym(2), that is, a 
Cartan subalgebra, consists of purely off diagonal matrices 


A = 


0 a 
a 0 


Oi G 


(4.3) 


A complementary subspace is Sym^(2), the space of diagonal 2x2 matrices. Notice 
that for any X £ Sym(2) we have 


[A,X]n = 


—2aa 0 
0 2ad 


(4.4) 


and hence, also in accordance with general theory, if a 7 ^ 0, then ad^ : Sym^(2) —> 
Sym^(2) is an isomorphism. Thus the inverse ad'^^ : Sym^(2) —> Sym^(2) is defined 
and hence 


ad^^ {ads X) = 


/9 


a 


a 0 
0 d 


for A = 


0 a 
a 0 


B = 


0 (3 

P 0 


, a AO. (4.5) 
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An operator of this form is called a sectional operator in the sense of Mishchenko 
and Fomenko [1976]. The equations defined by a sectional operator are 


A:= [A,ad^' (adBX)]^ = ^ 


—2ah 0 
0 2bd 


We shall now prove that (1.1) is not in this family. Indeed, since 


0 a + d 
—a — d 0 


XN + NX = 
equation (1.2) becomes 

A = [A, XN + AA] = {a+ d) 


= {a + d)N 


—2b a — d 
a — d 2b 


(4.6) 


(4.7) 


(4.8) 


The only way equations (4.6) and (4.8) can be identical is if one requires that a = d, 
which is not allowed since A is arbitrary in Sym(2). Therefore the system (1.1) is 
not in the list of equations of generalized rigid body type on 5p(2,M) described by a 
sectional operator in Mishchenko and Fomenko [1976]. 

Despite the fact that our system is not in the class of integrable systems studied 
in Mishchenko and Fomenko [1976], we shall see in the next sections that by using the 
techniques of Manakov [1976] and Magri [1978] (the method of recursion operators), 
the system is nonetheless integrable. 


5 Lax Pairs with Parameter 

To prove that system (1.1) is integrable for any choice of N, we will compute its flow 
invariants. Bear it in mind that, by virtue of the isospectral representation (1.2), 
we already know that the eigenvalues of A, or alternatively, the quantities trace A^ 
for A; = 1, 2,..., n — 1, are invariants. 

One way to compute additional invariants is to rewrite the system as a Lax pair 
with a parameter. One can do this in a fashion similar to that for the generalized 
rigid body equations (see Manakov [1976]). 

Theorem 5.1. Let X be a real parameter. The system (1.2) is equivalent to the 
following Lax pair system 

^{X + XN) = [X + XN,NX + XN + XN^] (5.1) 

Proof. The proof is a computation. The only nontrivial power of A to check is the 
first one. In fact, the coefficient of A on the right hand side of equation (5.1) is 

[A, NX + AA] + [A, A^] 

= A^A + AAA - AAA - AA^ + AA^ - A^A = 0, 


which proves (5.1). 
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We recall from Manakov [1976] and Ratiu [1980] that the left-invariant general¬ 
ized rigid body equations on SO{n) may be written as 

M = [M, 0], M(0) = Mo e 50(n), (5.2) 

where Q = Q~^Q € so(n) is the body angular velocity, Q G SO{n) denotes the 
configuration space variable (the attitude of the body), and 

M = J(0) := An -h OA G so(n) 

is the body angular momentum. Here J : 5o(n) ^ 5o(n) is the symmetric, positive 
definite (and hence invertible) operator defined by 

j(n) := An nA, 

where A is a diagonal matrix satisfying A* -|- Aj >0 for all i ^ j. For n = 3 the 
elements of Aj are related to the standard diagonal moment of inertia tensor I by 
/i = A 2 -I- As, /2 = As -I- Ai, Is = Ai -I- A 2 . 

Manakov [1976] has noticed that the generalized rigid body equations (5.2) can 
be written as a Lax equation with a parameter in the form 

^(Af+ AA2) = [M + AA2,n +AA]. (5.3) 

Note the following contrast with our setting: in the Manakov case the system 
matrix M is in 50 (re) and the parameter A is a symmetric matrix while in our case 
X is symmetric and the parameter N G 5o(re). 

For the generalized rigid body the nontrivial coefficients of A*, 0 < i < A: in the 
traces of the powers of M -|- AA^ then yield the right number of independent integrals 
in involution to prove integrability of the flow on a generic adjoint orbit of SO{n) 
(identified with the corresponding coadjoint orbit). The case i = 0 needs to be 
eliminated, because these are Casimir functions. 

Similarly, in our case, the nontrivial coefficients of A*,0 < i < k, in 

/i^(A:) ;= itrace(X +AiV)^, A: = 1, 2,..., re - 1 (5.4) 

yield the conserved quantities. The coefficient of A'',0 < r < k, in (5.4) is 

trace ^ r = 0,...,k, k = l,...,n-l, 

\i\=k-r\j\=r 

where i = ■ ■ ■ is), j = {ji, j 2 , ■ ■ ■ js) are multi-indices, iq,jq = 0,1,... ,k, and 

|A| = '^q^iiq, Ij] = Ylq=ijq- ^he coefficient of A^ is the constant so it should 
not be counted. Thus we have r < k. In addition, since the trace of a matrix equals 
the trace of its transpose, X G Sym(re), and N G 5o(re), it follows that 

trace = (-1)1^1 trace • • • X^^N^^X^K 
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Therefore, if r is odd, then necessarily 


trace 




^ = 0 


\i\=k-r\j\=r 


and only for even r we get an invariant. Thus, we are left with the invariants 

hk, 2 r{X) := trace E E (5.5) 

|'i|=A:—2r |j|=2r 


for k = 1,..., n — 1, iq = 1,..., fc, jq = 0,... ,k — 1, r = 0,..., , where [1] 

denotes the integer part of £ G M. 

The integrals (5.5) are thus the coefficients of A^'’, 0 < 2r < /c, in the expansion 
of ^trace(X + XN)^. For example, if /c = 1 or A: = 2 then we have one integral, 
the ceofficient of A^. If A: = 3 or A: = 4, only the coefficients of A^ and A^ yield non¬ 
trivial integrals. If A: = 5 or A: = 6 it is the coefficients of A^, A^, and A^ that give 
non-trivial integrals. In general, for the power k, we have integrals. Recall 

that A: = 1,..., n — 1. If n — I = 2.^, we have hence 


I + 1 + 2 + 2 + --- + 


n — 1 -|- I 


+ 


n — I -|- 1 


= l + l + 2 + 2 + ---+£ + i 


n — 1 I n — 1 


+ 1 = 


= £{£ + !) = 

integrals. If n — 1 = 2.^ -|- 1 then we have 

n — 2 -|- 1 

1 + 1 + 2 + 2 + --- + 


n — 1 n -|- 1 


+ 


ra — 2 -|- 1 


+ 


n — 1 -|- 1 


= I + I + 2 + 2 + ---+.f + ^ + (^ + l) 

= £{£+l) + {£+l) = {£+l)^=(!iy 

integrals. However, 


re 


re -1- 1 

2 


2 


' n — 1 re -|- 1 



if re is odd 
if re is even 


Concluding we have 


re 


re -1- I 

2 


2 


invariants which are the coefficients of A^^, 0 < 2r < A;, in the expansion of 
^ trace(X -|- XN)^ for A: = 1,..., re — 1. 

Are these integrals the right candidates to prove complete integrability of the 
system X = A^]? 
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If N is invertible, then n = 2p and hence 


n 

n + 1 


2p 

'2p + 1' 

2 

2 


2 

2 


1 

2 


+ p-p) 


= i (dimsp(2p,] 


ranksp(2p, M)) 


which is half the dimension of the generic adjoint orbit in 5p{2p, M). Therefore, 
these conserved quantities are the right candidates to prove that this system 
is integrable on the generic coadjoint orbit of Sym(n). This will be proved in 
the next sections. 


• If is non-invertible (which is equivalent to d ^ 0), then n = 2p -\- d and 
hence 


n 

n + 1 


2p + d 

2p -\- d+l 

2 

2 


2 

2 


= (P+ 2 

= p" +p 
= p^ +pd + 


p + 


d + 1 


+ 


d + 1 


+ 


'd' 

d+l 

2 

2 


d + l 
2 


The right number of integrals is p{p + d) according to Proposition 3.2, so this 
calculation seems to indicate that there are additional integrals. The situation 
is not so simple since there are redundancies due to the degeneracy of N. Note, 
however, that if d = 1, then we do get the right number of integrals. We shall 
return to the study of the degenerate case in §7. 


6 Involution 

In this section we prove involution of the integrals found in the previous section for 
arbitrary N G so(n). 

Bi-Hamiltonian structure. We begin with the following observation. 

Proposition 6.1. The system X = X‘^N - NX^ is Hamiltonian with respeet to the 
bracket {f,g}]^ defined in (3.1) using the Hamiltonian h 2 {X) := itrace(X^) and 
is also Hamiltonian with respect to the compatible bracket {/, g}p^ defined in (3.2) 
using the Hamiltonian hfiX) := ^ trace(X^). 

Proof. We have implicitly checked the first statement already using Euler-Poincare 
theory, but here is a direct verification. We want to show that the condition / = 
{/,/i 2 } 7 Y for any / determines the equations X = X'^N — NX‘^. First note that 
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/ = trace(V/(X)X). Second, since Vh2{X) = X, the right hand side {/,/i 2 }jv 
becomes by (3.1) 


{/, ^2}Ar {X) = - trace 
= — trace 


x(vf{X)NX - XNVf{X)^ 
(vf{X)NX‘^ - Xf{X)X‘^Ny 


Thus, X = X'^N — NX'^ as required. 

To show that the same system is Hamiltonian in the frozen structure, we proceed 
in a similar way. Noting that Xh^{X) = we get from (3.2) 

{/, {X) = - trace (vfNX^ - X^NXf) 

= - trace (vfNX^ - VfX^N ^, 


and hence X = X'^N - NX^, as before. 


Involution. Next we begin the proof that the [^] integrals given in (5.5), 

namely 

hk, 2 r{X) := trace E E 

\i\=k—2r \j\=2r 

where k = 1,... , n — 1, iq = 1,..., fc, jq = 0,... ,k — 1, r = 0,..., , are in 

involution. It will be convenient below to write the expansion of starting with 
the highest power of A, that is, 

k 

hl{X) = I trace (X + AN)'^ = V \^-'^hk,k-r{X). (6.1) 

Ki 

r=0 

As explained before, not all of these coefficients should be counted: roughly half of 
them vanish and the last one, namely, hk^ki is the constant . Consistent with our 
notation for the Hamiltonians, we set hk = hk^- 

Firstly we need the gradients of the functions h^. 

Lemma 6.2. The gradients V/i^ are given by 

Vhl{X) = ]^{X + AiV)^-i + i(A - ANf-^. (6.2) 

Proof. We have for any Y £ Sym(n), 

((V/i^(A),y)) = dhl{X) ■ Y = trace ((A + AXf-^Y'^ 

= i trace ( ((^ + AA)^"^ + (A - AA)^"^) p) . 


Since ((,)) is nondegenerate on Sym(n), the result follows. 
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Proposition 6.3. 

BxiVhUX)) = Cx{Vh^,^,{X)) (6.3) 

Proof. By (3.3) we have 

Bx{Xh^{X)) = XVhl{X)N - NVht{X)X 

= i X{X + XN)^-^N + X{X - XN)^-^N 

-N{X + AiV)^“^X - N{X - XN)^-^X 

= ^ {X + XN)’^N - XN{X + XN)^-^N + {X - XN)^N + XN{X - AiV)^“^iV 

-N{X + XN)^ + XN{X + XN)'^-^N - N{X - XN)^ - XN{X - XN)’^-^N 

= ^ {X + XN)^N + {X - XN)’^N - N{X + XN)’^ - N{X - XN)’^ 

= Xhl^,{X)N - iVVhti(^) = Cx{Xhl^,iX)) 

by (3.4), which proves the formula. ■ 

Proposition 6.4. The functions hk^k-r satisfy the recursion relation 

Bx{Vhk,k-r{X)) = Cx{Xhk+l,k-r{X)) (6.4) 

Proof. Substituting (6.1) into (6.3) we obtain 

k k-\-l 

X^-^Bx {Xhk,k-r{X)) = Y A"+'-''C'x {Vhk+l,k+l-r{X)) . 
r=0 r=0 

Since Xhk+i^k+iiX) = formula (3.4) implies that Cx {Xhk+i ,k+i{X)) = 0. 

Thus on the right hand side the snm begins at r = 1. Changing the summation 

index on the right hand side now from r to r — 1 and identifying the coefficients of 

like powers of A yields (6.4). ■ 

Remark. It is worth making a few remarks about Propositions 6.3 and 6.4. Note 
that unlike the similar recursion for the rigid body Mankov integrals (see e.g. Morosi 
and Pizzocchero [1996]) our polynomial recursion relation (6.3) does not have a 
premultiplier A on the right hand side and the polynomials on the left and right 
hand sides appear to be of different order. This cannot be and indeed is not so. 
Indeed the highest order order coefficient on the right hand side vanishes by virtue 
of following result. 

Corollary 6.5. The functions hk^k-i{X) are Casimirs for the frozen Poisson struc¬ 
ture, i.e. 

Cx{Xhk,k-i{X)) = 0 (6.5) 


for all k. 
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Proof. By (6.1), hk,k-i{X) = trace so its gradient equals V/ifc,fc-i(^) = 

So (3.4) immediately gives (6.5). ■ 

The recursion relations (6.4) for r = 0 also imply the following relation between 
the Hamiltonians that can also be easily checked by hand. 

Corollary 6.6. 

Bx {Vhk{X)) = Cx (Vhfc+i(X)) (6.6) 

Example: An interesting nontrivial example of the recursion relation to check 
is Bx{dh^^ 2 {X)) = Cx{dhi^ 2 {X)) where h^^ 2 {X) = trace(iV^A) and /i 4 , 2 (-A) = 
trace(A^^X^) + ^ trace(A^XiVA). This example illustrates how the recursion relation 
works despite the apparent inconsistency in order. 

Uising the recursion relations involution follows immediately. 

Proposition 6.7. The invariants hk^k-r in involution with respect to both Pois¬ 
son brackets {f,g}x and {f,g}px- 

Proof. The definition of the Poisson tensors Bx and Cx and the recursion relation 
(6.4) give 

{hk,k—rthl^l—q\j^ — {(S/ hk^k—r{X'), Bx {X 

= {{Xhk,k-r{X),Cx{Xhl+l,l.q{X)))) 

\hk,k—r) ^Z+1,Z— {dl+l,l—qi hk,k—r\px 
= -{{Vhl+l,l.q{X),Cx{Xhk,k-r{X)))) 

= -{{Vhl+l,l.q{X),Bx{Vhk-l,k-r{X)))) 

— \kil+l,l—q)hk—\^k—r\x — {^fc—1,A:—rj 

for any k,l = 1,...,re — 1, r = 1,... ,k and g = 0,— 1. Of course, in these 
relations we assume that k — r and I — q are even, for if at least one of them is odd, 
the identity above has zeros on both sides. Repeated application of this relation 
eventually leads to Hamiltonians hk^k-r where either k — r is a power of A that does 
not exist for k, in which case the Hamiltonian is zero, or one is led to /iq.o which is 
constant. This shows that {/ifc.fc-n ~ ^ for any pair of indices. 

In a similar way one shows that {hk,k-r, ~ * 


7 Independence 

To complete the proof of integrability we need to show that the integrals hk^ 2 r are 
independent. We will demonstrate this first in the generic case that N is invertible 
with distinct eigenvalues. 

By (5.5), the gradients of the integrals /ifc, 2 r have the expression 

V/lfc,2r(Ai) : = E E 

2r—1 |j|=2r 

where A; = 1,..., re - 1, ig = 1,..., A:, jg = 0,..., A: - 1, r = 0,..., [^]. 


(7.1) 
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The Generic Case. We consider the case N invertible with distinct eigenvalues. 
Therefore d = 0 and n = 2p. In this case we show that the integrals hk^ 2 r given in 
(5.5) are independent, and hence the system (1.1) is system is integrable. 

Theorem 7.1. For generic N the integrals hk^ 2 r given by equation (5.5) are inde¬ 
pendent. 

Proof. We are concerned with the linear independence (in a generic sense) of (7.1) 
where k = I,... ,n — 1, iq = 1,... ,k, jq = 0,... ,k — 1 and r = 0,... [^(/c — 1)]. We 
assume that N is invertible with distinct eigenvalues and, without loss of generality, 
that X is diagonal, 

X = diagp. 

This reduces the problem to a problem about the independence of polynomials 
in single matrix variable. 

Now, we aim to prove a stronger statement: the terms 

Vij = 

are independent for all multi-indices i and j in the above range. Note however 
that each Vij is a g-degree polynomial in pi, p 2 , ■ ■ ■, hn, where q = k — 2r — 

{0,..., n — 2}. Let 

T-Cg = {vij : |i| = q, |j|even}. 

Clearly, in a generic sense, if linear dependence exists, it must exist within a set 
7iq. In other words, if we can prove that there is no linear dependence within each 
7iq, we are done. (Note that since A: < n — 1 in the expression (7.1) there is no 
dependence of powers of X on lower powers through the characteristic polynomial.) 
There is nothing to prove for q = d For <7 = 1 we have 

Til = {XN^ : j even} U X : j even}. 

Suppose that there exists linear dependence in Tii. Then there necessarily exist 
Pq,P 2 ,..., Pn -2 and kq; ■ ■ ■ > not all zero, such that 

X (E (E X = 0 = XR{N) + K{N)X = 0. 

Therefore, 

Pa[R{N)]a,h + [K{N)]a^kPi, = 0, o, 6 = 1,..., n. 

Generically (i.e., for all p except for a set of measure zero) this can hold only if 
R{N), K{N) = 0. But deg 7?, deg X < n — 1 and, since the eigenvalues of N are 
distinct, the degree of the minimal polynomial of N is n. Therefore K,R = 0, a 
contradiction. Hence there is no linear dependence. 

We continue to s = 2. Now 

U 2 = {X*iX^'iX*2X^2X*3 : + Z 3 = 2 , ji + j 2 even}. 
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Assume that there exist not all zero, s.t. 

= 0 . 




Therefore 


^ Y, = 0, a, c = 1,..., n. 




Note that we want the above to hold for all real pk, but this is possible only if 






thus 


= 0 . 




We again obtain a polynomial in N'^ of degree < n/2, which cannot be zero: a 
contradiction. 

We can continue for higher s in an identical manner. ■ 

Hence, since we have involution and independence we have proved the following. 

Theorem 7.2. For N invertible with distinct eigenvalues the system (1.1) is com¬ 
pletely integrable. 

Corollary 7.3. For N odd with distinct eigenvalues and nullity one, the system 
(1.1) is completely integrable. 

Proof. In this case we have d = 1 and n = 2p + 1. All eigenvalues are distinct with 
one of them being zero. The above proof of indepdence still holds, the only change 
being that the characteristic (and mininal) polynomial of N is of form Nw{N‘^), 
where rc is a polynomiail of degree (n — l)/2. ■ 


8 Conclusion and Future work 

We have demonstrated integrability of the system (1.1) for appropriate N by showing 
involution and independence of a sufficient number of integrals. It is also of interest 
to analyze linearization on the Jacobi variety of the curve 

det(zl - AA - A) = 0 

using the theory discussed in Adler, van Moerbeke, and Vanhaecke [2004] and Grif¬ 
fiths [1985], for example. We shall discuss these algebro-geometric aspects in a 
future paper. Independently Li and Tomei [2006] have shown the integrablity of the 
same system in precisely the two cases discussed in this paper employing different 
techniques; they use the loop group approach suggested by the Lax equation with 
parameter (5.1) and give the solution in terms of factorization and the Riemann- 
Hilbert problem. 

Another interesting variation of this system that we shall consider in future work 
is the following. 
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A generalized system. The flow of (1.1) can be rendered more general by 
complexiflcation. Generalizing it to evolution in su(n) yields an n^-dimensional 
flow of generalized rigid body type with two natural Hamiltonian structures. Let 
Xq G su(n), N G Sym(n,]R), and consider 

X = [X^,N] = [X,XN + NX], X{0)=Xo. (8.1) 

Note that X(t) evolves in su(n) since one readily checks that [X, XN+NX] G su(n). 
Moreover, one can generalize this still further and take N G su(n). We define 

Hi{X) = ^ traceX{XN + NX), 

H 2 {X) = ^ trace 

Note that both Hamiltonians are real and that H 2 gives us our earlier Hamiltonian 
in the case that X is symmetric but that Hi is zero in this case. 
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